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ABSTRACT 

Ultraluminous X-ray sources (ULXs) have been described variously as the most luminous normal 
X-ray binaries, as hypernovae, and as "intermediate-mass" black holes with masses of hundreds to 
thousands of solar masses. We present results on M74 X-l (CXOU J013651. 1+154547), a ULX in the 
nearby spiral galaxy M74 (NGC 628), from observations by Chandra and XMM-Newton. M74 X-l is 
a persistent source which exhibits extreme variability and changes in spectral state on timescales of 
several thousand seconds. Its variability timescales and behavior resemble some Galactic microquasars. 
This suggests that the emission mechanism may be related to relativistically beamed jets, and that 
M74 X-l could be an extragalactic "microblazar" — a microquasar whose jet axis is aligned with our 
line of sight. We also note that its spectrum is consistent with the presence of a low-temperature disk 
blackbody component, which, assuming it is due to radiation from an accretion disk, could indicate 
that M74 X-l contains an intermediate-mass black hole. 

Subject headings: galaxies: individual (NGC 628, M74) X-rays: individual (CXOU 

J013651. 1+154547, M74 X-l) 
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1. INTRODUCTION 

Ultraluminous X-ray sources (ULXs) are historically 
defined as point sources whose luminosities exceed the 
Eddington luminosity for a 10 M Q compact object, L x ~ 
10 39 erg s -1 . Although th ey were discovered in the 
Einstein era i|Fabbiandll989|) . and more than 100 ULX 
candi dates were subsequently ca t alogued by the ROSA T 
HRI l|Roherts fc Warwick! 1200(1 IColhert fc Ptakl 12002ft . 
much about the nature of these sources was unknown 
until recent times. Observatories such as ASCA pro- 
vided val uable information abou t the spectra of ULXs 
(see, e.g.JMakishi ma et alJ 2000): however, except in the 
cases where the sources were significantly variable, there 
remained the possibility that what appeared to be a ULX 
was really a group of spatially unresolved sources. 

The ability to spatially separate ULXs from their sur- 
rounding environs is now possible with the superb res- 
olution of the Chandra X-ray Observatory, and both 
Chandra and XMM-Newton can perform high quality X- 
ray spectroscopy and temporal analysis. This affords 
us detailed insight into the nature of ULXs: we now 
know, for example, that probably 10% or fewer ULXs 
are a ssociated with recent supernovae ijRoberts et al.l 
2002), and that the others are, in fact, isolated sources. 
Although some ULXs are found in elliptical galaxies 
l)Colbert fc Pta k 2002), most reside in star forming and 
starburst galaxies, implying that they are associated with 
recent star formation (though the exact nature of this 
link is unclear). 

A number of models have been proposed to explain 
the high X-ray luminosity of ULXs, and it is likely 
that more than one is required to explain all types 
of observed ULX acti vity. Some ULXs may be extra- 
galactic microquasars ijkordin g, Falcke fc Markoflj|200 
|Fojichnri_^^ Aharonian fc Kir' 

2002; King et al. 2001), in which case the observed high 
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luminosities would be due to relativistic beaming. Oth- 
ers may be stellar-mass blac k holes with super-Eddington 
radiation ( Be gelmanl 1200 21 . Some are almost certainly 
background AGN that happ en to fall along the line- 
of-sig ht to the host galaxy l|Irwin. Bregman. fc Athevl 
2004). One of the most intriguing suggestions is that 
t he y are intermediate-mass black holes (IMBHs, see 
e.g. IColbert fc Mushot^klll999l: iMakishima et al.ll2000t 
Mille^SHl003). Aside from their high luminosities, 
the low-temperature blackbody component in some ULX 
spectra can be used to argue for the presence of an IMBH 
( Mil ler et alJ2003|) . In this scenario, the blackbody spec- 
trum is assumed to arise from an accretion disk, and its 
temperature can be used to infer the radius of the disk's 
inner edge and thus a mass for the central object. 

In this paper we focus on observations of a specific 
ULX, M74 X-l (CXOU J013651. 1+154547). The first X- 
ray imagi ng observation of M74, which is at a distance of 
8.7 Mpc ijHuchra. Vogelev fc GellerHTflfl^ . was obtained 
with the Einstein HRI on 1980 January 5. Though it 
is not included in the Einstein source catalog, there is a 
source present at the location of M74 X-l with a 0.3-10 
keV flux of 7.6 x 10~ 14 erg cm" 2 s" 1 (14 total counts), 
corresponding to a luminosity of 7 x 10 38 erg s . M74 
was also observed by ASCA on 1995 July 28, but no 
source is detected at the location of M74 X-l. We 
derive an upper limit to the 0.2-10 keV flux of 7 x 
10" 14 erg cm" 2 s _1 (corresponding to a luminosity of 
6 x 10 38 erg cm" 2 s" 1 ). No pointed ROSAT observa- 
tions were performed. 

The Chandra observations of M74 were conducted as 
part of a Cycle 2 large project to study the p oint source 
popu lation of nearby, face-on spiral galaxies ijPrestwichl 
200 lj. M74 X-l, one of the objects in a subsample of 
ultraluminous point sources, was discovered to be ex- 
tremely variable and exhibit strong flares with a char- 
acteristic timescale of several thousand seconds. In Sec- 
tion we report on the recent X-ray observations of 
M74 X-l; in particular, its spectral characteristics and 
variability. In Section |21 we present optical and radio ob- 
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servations of M74. In Section 0] we discuss M74 X-l's 
location in M74, as well as its observed spectral and vari- 
ability properties; and compare our findings with Galac- 
tic microquasars and with other ULXs. 

2. X-RAY OBSERVATIONS AND DATA ANALYSIS 

Chandra has observed M74 twice, on 2001 June 19 and 

2001 October 19. Each observation utilized the ACIS 
back illuminated S3 chip, and basic processing was done 
by the Chandra X-ray Center (CXC). There are two small 
background flares in the first observation during which 
background rates jump to two and three times the quies- 
cent level. These flares were not filtered out for analysis 
because they contribute a negligible (~ 0.3) number of 
counts to the source data presented here. The total du- 
ration of the "good time intervals" (GTIs) for the June 
dataset is 46.4 ks. There are no background flares in the 
October observation, whose GTIs sum to 46.2 ks. Both 
datasets have been filtered to include only energies in the 
range 0.3-10.0 keV. 

XMM-Newton observed M74 for 36 ks on 2002 Febru- 
ary 2. The data from MOS1, MOS2, and the PN were 
reduced individually, then fit simultaneously during spec- 
tral analysis, but were summed together for the temporal 
analysis. Data reduction was performed using SAS (Sci- 
ence Analysis Software) v5.3. Several significant back- 
ground flares were removed for analysis, resulting in the 
loss of ~ 7 ks of data. The data were also filtered for 
"good" events plus pattern < 12 (MOS) and < 4 (PN), 
and for the energy range 0.2-10.0 keV. The total dura- 
tion of the GTIs are, for MOS1 and MOS2, 28.9 ks; and 
for the PN, 28.7 ks. 

2.1. Spectral Fitting and Energetics 

Spectral fits were performed using Sherpa, the fitting 
tool included in the CIAO data analysis package. Back- 
ground data were subtracted prior to fitting, powell 
was used for optimization, and Chi Data Variance was 
the statistic function employed. The fits to the Chan- 
dra (XMM) data were performed on data grouped in 15 
(20) count bins. The Chandra ARFs were corrected for 
the ACIS low energy QE degradation using the corrarf 
program. 4 In the XMM observation, M74 X-l is 6' off- 
axis and close to chip gaps on all the detectors, so the 
tools rmf gen and arf gen were used to produce correct 
response files for fitting. 

For the 2001 June 19 (Chandra) observation, there are 
a total of 347 source counts, of which approximately 
seven are due to background contamination. The 2001 
October 19 (Chandra) observation contains a total of 963 
source counts, of which approximately eight are due to 
background. Both source ellipses are of a size that is ex- 
pected to contain 98% of the source energy, given their 
off-axis angle (3.'6 for June, and 3f 8 for October). For the 

2002 February 2 (XMM-Newton) observation, a source- 
free background region on the same chip as the source 
was used for background determination. After back- 
ground subtraction, there are a total of ~ 650 counts in 
the PN spectrum and ~ 430 counts in the MOS spectra. 

At its peak brightness in the Chandra observations, the 
0.3-10.0 keV count rate is 0.11 cts s" 1 , corresponding to 

4 Information about corrarf can be found at 
[http://cxc.harvaxd.edu/cal/Links/ Acis/acis/Cal-prods/qeDeg/| 
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Fig. 1. — 2001 June 19 Chandra spectrum with the absorbed 
(Galactic Nn) power-law fit overplotted and residuals shown below. 



a source flux of 1.3 x 10~ 12 erg cm~ 2 s -1 , or a luminos- 
ity of 1.2 x 10 40 erg s _1 . To verify that the spectrum is 
not affected by pileup, we used the CIAO tool mkpsf to 
create a PSF image at the location of M74 X-l. The inner 
3x3 pixel region of the PSF contains 43.4% of the total 
flux, whereas this would be 88.6% were the source at the 
aimpoint of the detector. Using the pileup algorithm im- 
plemented in PIMMS (Portable, Interactive Multi-Mission 
Simulator), the source at its brightest (0.11 cts s _1 ) is 
3.2% piled, while on average (a rate of 0.04 cts s~ 4 ) it is 
1.2% piled. This will not cause a statistically significant 
change in the shape of the spectrum. 

For initial spectral fitting, we tried a range of standard 
simple models, such as a disk blackbody, power-law, ther- 
mal bremsstrahhmg, and MEKAL. Of these models, only 
the power-law provides acceptable fits (see Figures Q [21 
and EJ. However, the 2001 October 19 spectrum still 
has significant residuals beyond the power-law fit. If 
we restrict the fitting range of the power-law to cover 
1.5-10.0 keV, these residuals are manifested as a soft 
excess (see Figure |2J. We were motivated by the spec- 
tra se^njii_Galac_tic_X-ray binaries (and in other ULXs, 
c.f. IMiller et all 12003]) to add a disk blackbody (DBB) 
component to the power-law, and using an F"-test, we 
find that this improves the fit at the 3.3<r level (see Fig- 
ure [Sj . We find that the soft excess is similarly well- 
fit by a MEKAL model with low abundances and by 
a single-temperature blackbody spectrum. Using an F- 
test to compare these models with the power-law fit, we 
determine that the addition of a MEKAL component is 
favored at the 3.7<7 level, whereas the blackbody compo- 
nent is favored at the 3.1cr level. We also fit the spec- 
trum with a broken power-law model, but find that the 
parameters are not well constrained. We present the fit 
parameters for all the datasets in Tables ^ (power-law 
and power-law plus disk blackbody fits) and [2] (power- 
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TABLE 1 

Power-law and Power-law Plus Disk Blackbody Spectral FiTS a 



Observation 






Power-law 




Disk Blackbody 




xl (DOF, sig.) 




(xio 2D ) 


r 


Amplitude 13 

(xio D ) 


Flux c 

(xlO ld J 


(keV) 


Norm 


Flux c 

(xio ±J ) 




2001 June 19 


4.8 


9 9C -+0.19 
z - zo -0.18 


0.96 ±0.09 


0.49 








1.67 (19, N/A) 




4±i 3 

4.8 

7+16 
1 -7 


7 +0.9 
u -'-1.6 


O qo+0-52 
u - yz -0.20 

17+°- 3 l 
u - 1 '-0.14 
n iq+0-38 
u - iy -0.17 


0.48 
0.45 
0.45 


0.29+°-{£ 

— U.UO 

27+ - 13 
u - z '-0.10 


0.22±° H 

u - c>z -0.26 


0.25 
0.27 


1.75 (18, 0.2cr§ ) 
1.25 (17, 2.1ct§ ) 
1.32 (16, 2.0(7+ ) 


2001 Oct 19: All 


4.8 
4.8 

^.4-12 
6 


2.02 ±0.10 

, OO+0.21 
1-»»_0.17 

1 I" 1 " 04 
1 - i -0.5 
^ -4-0 48 

l- lo -0.71 


2.72 ±0.15 

2 45+0- 38 
z '^°-0.28 

1 0+°' 6 

-i nr,4-0 87 


1.52 
1.49 
1.35 
1.36 


0.27 ±0.05 
0.26±g;lg 


6 +0 - 6 
u '"-0.3 

0.8lJ 2 6 8 


0.54 
0.56 


1.23 (55, N/A) 
1.22 (54, l.lcr§ ) 
0.96 (53, 3.5o- § ) 
0.98 (52, 3.3o-t ) 


Bright: 


4.8 
q+2.8 
U -JV/A 

4.8 

-,,+19 
±J -13 

4.8 


, 44I 17 

i-58lS:iI 

67+ ' 50 
u -°' _o 61 

n oc+0.54 

0.83«;f 8 t 


1 O+0.47 

5-18_q 4g 

4 44+0-54 

1 6+ 14 
r- u _0 9 

2 +2 '° 


3.39 
3.44 
3.97 
3.93 
3.65 


29+ 08 
u - za -0.07 

22+ ' 18 
u - z -0.08 

n si+o- (l H 

U - ,}i -0.05 


1 0+ 18 

4 q+390 

fiq+0'58 
U ' DJ -0.30 


1.09 
1.7 
1.01 


1.70 (27, N/A) 
1.58 (26, 1.7cr§ ) 
1.16 (25, 2.9cr§ ) 
1.19 (24, 2.5(7+ ) 
1.20 (25) 


Faint: d 


4.8 

7 i +8.9 
' -1-7.1 

4.8 


2.29i{j;£ 
0.83±°;i|t 


1 co+0.14 

0.29±° : ?? 


0.93 
0.99 
0.58 


n 0,4-0.05} 
u - ,31 -0.05 


n OO+0.34 
u -°°-0.16 


0.48 


0.86 (26, N/A) 
0.88 (25, 0.6cr § ) 
1.13 (26) 


2002 Feb 2 


4.8 


l- au -0.09 


2.38 ±0.30 


1.4 








1.09 (59, N/A) 




1.1 ±2.7 


1.77 ±0.17 


2.09 ± 0.34 


1.1 








1.04 (58, 1.9cr§ ) 




4.8 


1.35 ±0.39 


1.3 ±0.7 


1.3 


0.25 ±0.06 


0.6 ±0.6 


0.3 


0.92 (55, 2.6cr§ ) 




1.2 ±6.3 


1.24 ±0.52 


1.1 ±0.8 


1.3 


0.31 ±0.16 


0.2 ±0.4 


0.3 


0.92 (54, 2.1o- 1 ) 



a Upper and lower limits are the 90% confidence range. TVh was frozen to the Galactic value of 4.8 X 10 20 cm 2 for fits where there is no TVh 
confidence range given; otherwise, JVh was allowed to vary without bound. 
^Thc amplitudes of the power-laws are the flux at 1.0 keV in units of photons keV -1 cm -2 s — 1 , 
c 0. 3-10.0 keV flux in units of erg/cm 2 /s. 

^Adding a disk blackbody component actually increased the value of xt f° r these data. 

^Significance is determined using an i^-test comparing this fit to a power-law fit of the same data with TVh allowed to vary. 
^Significance is determined using an F-test comparing this fit to a power-law fit of the same data with Nu frozen to the Galactic value. 
''"For these fits, the "faint" and "bright" data were fit using a common power-law slope and blackbody temperature, allowing only their 
normalizations to vary. 
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Fig. 2. — 2001 October 19 Chandra spectrum with the absorbed 
(Galactic A^h) power-law fit overplotted and residuals shown below. 



law plus MEKAL fits). 

To explore possible spectral variability in the tempo- 
rally variable October data (the lightcurve is presented 
in Figure HJ), we separated it into "bright" (count rate 
> 0.03 cts s _1 ) and "faint" (count rate < 0.03 cts s _1 ) 
segments. Both segments had comparable numbers of 
counts; 490 in the bright data and 469 in the faint data, 
allowing for similar statistical properties. (There are an 
insufficient number of counts in the June observation 
to allow for a similar separation.) Fitting with simple 
power-laws, we find that the bright data are somewhat 
harder than the faint data: the best-fit photon index for 
the bright data is 1.6^5" l5 whereas for the faint data it is 

2.4^Q4' We also note that there appears to be a soft ex- 
cess present in the bright data — as with the full dataset, 
this can be modeled variously with a disk blackbody, 
MEKAL component, and single-temperature blackbody. 
Although the spectra appear to be different, it is pos- 
sible that the physical parameters of the system (for 
example, the power-law slope and blackbody tempera- 
ture) remain the same as the source brightens, and that 
only the contributions of different components vary. To 
test this, we fit the bright and faint data simultaneously 
with an absorbed power-law plus a disk blackbody. We 
used a common power-law slope and disk blackbody tern- 
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TABLE 2 

Power-law Plus MEKAL Spectral Fits 11 



Observation jVg Power-law MEKAL \l (DOF, sig.) 

T Amplitude 10 Flux c Plasma Temp. Abund. Norm Flux c 

(xio 20 ) (xio- 5 ) (xio- 13 ) (keV) (xio-6) (xio- 13 ) 



2001 June 19 


4.8 


i qq+0-33 
1,! * 5 — 0.3B 


0.76 ± 0.18 


0.45 


25+ 11 


i 


3.0+H 


0.06 


1.60 (17, 1.1<t§ ) 


2001 Oct 19: All 


4.8 
q+2.1 


1.93^ 

i.68±S:i; 


2 42+ u - M 
z - 4z -0.25 
1 QK+0.26 

Lyj-o.is 


1.43 
1.38 


7 +ii ' 2 

u -'-0.4 

7+ - 2 
u -'-0.2 


l 
l 


2.8ti; 8 7 
3.3 ± 1.5 


0.08 
0.09 


1.15 (53, 1.9cr§ ) 
1.03 (52, 2.9<rt ) 




4.8 
2.8+1* 


7 +0.5 

u -'-0.6 

0.7 ±0.6 


=+0-6 
U '°-0.3 

6+°' 7 
u ' D -0.4 


1.30 
1.29 


0.8 ±0.2 
0.8 ±0.2 


02+ 02 
u - uz -o.oi 

02+ - 04 


153+32 

126+ 95 


0.72 
0.64 


0.90 (52, 3.6<r!i ) 
0.90 (51, 3.7<rt ) 


Bright: 


4.8 


1 40+ U - 28 
1, * u — 0.29 


• 3 - D ' -0.87 


3.44 


0.7 ±0.2 


1 


11.7 ±5.7 


0.32 


1.40 (25, 2.1a* ) 




0+2.4 

"-N/A 


1.CSU_ 25 


3.31 ± 0.69 


3.47 


0.7 ±0.2 


1 


9.9l 4 4 9 7 


0.27 


1.24 (24, 2.3crt ) 




4.8 

0+8.1 


44+0-51 
u '^-0.63 

o w+o- 48 

U - OD -0.30 


1 1+ 1 ' 1 

i ' i -0.7 

1 3+ 1 - 1 


4.02 
3.93 


0.8 ±0.2 
0.8t°;» 


04+ 04 
06+ 06 


249+gg 
160t 2 °° 


1.29 
0.95 


0.98 (24, 3.3<r§ ) 
1.00 (23, 3.0crt ) 



2 for fits where there is no Nu confidence range 



a Uppcr and lower limits arc the 90% confidence range. TVh was frozen to the Galactic value of 4.8 X 10 2 
given; otherwise, Nu was allowed to vary without bound. 
^The amplitudes of the power-laws arc the flux at 1.0 keV in units of photons keV - 1 cm -2 s —1 . 

^Significance is determined using an .F-test comparing this fit to a power-law fit of the same data with TVh allowed to vary. 

§ Significance is determined using an F-test comparing this fit to a power-law fit of the same data with iVn frozen to the Galactic value. 
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Fig. 3. — 2002 Feb 2 XMM-Newton spectrum with the power- 
law model (Galactic iVjj) overplotted and residuals shown below. 
PN data are plotted with circles, MOS1 data are plotted without 
a symbol, and MOS2 data arc plotted with squares. 
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Fig. 4. — 2001 October 19 Chandra spectrum with a power- 
law fit over the restricted energy range 1.5-10 keV overplotted and 
residuals shown below. 



perature, and allowed only the normalizations to vary. 
Both datasets are consistent with a power-law slope of 
r = 0.8 ± 0.6 and a disk temperature of 0.3llg;°5 keV, 
but we must allow the amplitudes of both components 
to vary in order to get reasonable fits for the bright data. 
Plots of the "faint" and "bright" spectra are presented 
in Figures [7| and [SI and the results of spectral fitting are 
included in Tables and [21 

To further explore the spectral variability, we plotted 
lightcurves of M74 X-l in hard (1.0-10.0 keV) and soft 
(0.3-1.0 keV) bands. These lightcurves, along with a plot 



of the hardness ratio, are presented in Figure[51 Unfortu- 
nately, the limited number of counts (and therefore the 
statistical noise) hinders detailed interpretation. There- 
fore, we also binned the data according to the source 
count rate, and calculated its hardness ratio. There 
appears to be a trend toward increasing hardness with 
higher count rate, as can be seen in Figure ^| although 
a linear fit with non-zero slope is preferred over a flat 
distribution at only the 1.6a level. 

We performed a similar division of the 2002 February 2 
dataset, but find that spectral hardening with increasing 
count rate is not present in these data. 
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Fig. 5. — 2001 October 19 Chandra spectrum fit by a power- 
law plus disk blackbody. The fit and individual componenets are 
overplotted and residuals are shown below. 
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Fig. 6. — Lightcurves for all three observations. The data are 
binned by 500 s, and the full energy range (0.3-10.0 keV) has been 
used. 
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Fig. 7. — 2001 October 19 "faint" Chandra spectrum with the 
power-law model overplotted and residuals shown below. 
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Fig. 8. — 2001 Octoberl9 "bright" Chandra spectrum with the 
power-law model overplotted and residuals shown below. 



2.2. Variability 

M74 X-l is active and extremely variable in all three 
observations, with count rates sometimes increasing by 
an order of magnitude on a timescale of several thousand 
seconds. This variability is not caused or influenced by 
spacecraft dither or instrumental effects, and appears as 
a sequence of "flares" . The rise and decay of these flares 
generally appear symmetric, with a few exceptions when 
a flare is followed by a period of "noisy" decay (for ex- 
ample, after the fourth and seventh flares in the 2001 



October 19 data; see the lightcurves presented in Fig- 
ure EJ. The source is relatively inactive during the first 
half of the 2001 June 19 observation, but exhibits several 
significant flares in the latter half. It is active throughout 
the 2001 October 19 and 2002 February 2 observations; 
the former contains one particularly energetic flare dur- 
ing which the count rate rises from (0.004 ± 0.003) cts 
s" 1 to (0.08 ± 0.01) cts s" 1 in 2500 s. 

An FFT reveals that there are a number of signifi- 
cant extremely low-frequency peaks (which could be in- 
terpreted as quasi-periodic oscillations, or QPOs) in the 
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Fig. 9.— 2001 October 19 lightcurves in hard (1.0-10.0 keV) and 
soft (0.3-1.0 keV) bands, with the hardness ratio plotted below. 
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Fig. 10.— Hardness ratio (1.0-10.0 kev)/(0.3-1.0 keV) plotted 
against count rate for the 2001 October 19 data. Counts have been 
summed over specific count rate intervals in order to improve the 
statistics of the hardness ratios. A linear fit to the data shows that 
a line of slope 11 is preferred over a line of slope zero at the 1.7a 
level. 
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Fig. 11. — Normalized power density spectrum of the 2001 Oc- 
tober 19 observation. The FFT was performed with a bin time 
of 3.241 s, and the frequency bins were averaged to improve the 
statistics. 
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Fig. 12. — Period folding for all three observations. Folding was 
performed with the IRAF PROS routine period from 3000-9000 s 
with 50 s binning. 



power spectrum of the October observation (see Fig- 
ure lll|) . However, there are no significant peaks at higher 
frequencies up to the Nyquist frequency of 0.156 Hz in 
either of the Chandra observations, both of which have 
frame times of 3.2 s. Since the XMMPN data have time 
resolution of 73.4 ms (and a total of 660 counts), we 
searched it for higher frequency signals, but did not find 
anything significant. Performing period folding using the 
IRAF task period on all three of the observations pro- 
duces x 2 values for a range of periods (see Figure [TJJ) . 
There are broad peaks present in each, although there is 
not a single fundamental frequency. 

3. OPTICAL AND RADIO DATA 



M74 was observed by the Gemini North telescope as 
part of its commissioning campaign with the GMOS de- 
tector on 2001 August 13-14. Four images in each filter — 
G (475 nm), R (630 ran), and / (780 nm) — were obtained 
and co-added for a total of 960, 480, and 480 s respec- 
tively. Basic processing was done at the Gemini Obser- 
vatory. To register the world coordinate system (WCS) 
of the i?-band image, we first extracted bright stars us- 
ing the SExtractor program, then used the imwes rou- 
tine (both from the WCSTools package 5 ) in conjunction 
with the USNO-A2.0 catalog to match 93 M74 sources 
to USNO-A2.0 objects. Sources in the G- and /-bands 

5 WCSTools documentation can be found at 
|http: / / tdc-www~h arvard.edu /sof tware/wcstools / 1 
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Fig. 13. — Gemini GMOS data (GRI image courtesy Gemini Ob- 
servatory/GMOS team) overlaid with regions for OB associations 
118, 121, and 122 Hvanov et aljlT 9921 and the 0.2" positional un- 
certainty of M74 X-l (small circle). 



were then centroided and these images were aligned to 
the i?-band image using IRAF. After this process was 
complete, the centroids of three bright X-ray sources 
in the field which have optical counterparts fell within 
± 0'.'2 of them. Within the 0'.'2 Chandra error circle 
there is no counterpart down to the limiting magnitude of 
my — 25 (calculated using the GMOS North Integration 
Time Calculator 6 , assuming a detection signal-to- noise 
threshold of 5.) 

We also obtained observations of M74 taken on 2001 
December 16 with the Fred Lawrence Whipple Obser- 
vatory (FLWO) 1.2 m telescope. The galaxy was im- 
aged in UBVRI (Harris broadband filter set), SII, and 
Ha. We note that there is no obvious ionization nebula 
present in either the SII or Ha i mages, as is sometimes 
observed in the vicin ity of ULXs l|Pakull fc Mirion i 2002; 
iRoberts et alJl2003|) . A Very Large Array (VLA) obser- 
vation at 1.425 GHz did not reveal a ny radio counter - 
part with an upper limit of 80 /xJy l)Stockdalel 12002). 
The X-ray flux density of the 2001 October 19 spec- 
trum is 4 x 10~ 14 erg cm~ 2 s _1 kcV -1 at 1.0 keV, which 
gives a flux density of 1.7 x 10 -2 /iJy. Given this X-ray 
flux, the lower limit to the X-ray-to-radio flux ratio is 
log[F x /F R ] > -3.8. 

4. DISCUSSION 

M74 X-l has some very unusual observational char- 
acteristics. It is a ULX (peak luminosity of 1.2 x 

6 This calculator can be found at 



http: / / www. gemini . edu /sciops /instruments /gmos /gmoslndex. html 
following the .Performance and Use link. 



10 40 erg s _1 ) which is located near a young star clus- 
ter in M74. Its spectrum is reasonably fit by a simple 
power law, but in some observations there appears to be 
an excess at lower energies. This excess can be fit by a 
low temperature disk blackbody (0.3 keV) or a MEKAL 
model with low abundances. Its flux varies by an order 
of magnitude on timescales of half an hour, and there is 
tantalizing evidence for a QPO. The spectrum appears 
to harden as the source brightens, although the signifi- 
cance of this effect is low. In this section, we consider 
the physical interpretations for these observations. 

4.1. Location of M74 X-l 

A major consideration in the study of ULXs is whether 
or not they actually reside in the galaxy with which they 
are coincident or whether they are, in fact, background or 
foreground sources. If M74 X-l were in the Milky Way, 
given its Galactic latitude of —45.7 and assuming it is 1.4 
kpc above the Galactic plane (in the thick disk popula- 
tion), it would be at a distance of ~ 2 kpc. If it were an 
LMX B, its counterpart would like ly have 4 > My > — 4 
ijvan Paradiis fc McClintocklll994() . which at a distance 
of 2 kpc (and unobscured) would have 15.5 > my > 7.5. 
Although the optical and X-ray observations were not 
simultaneous, the source seems to be X-ray active from 
June 2001 to February 2002, and two of the optical ob- 
servations (Gemini and FLWO) were within this period. 
No object brighter than the my = 20.6 OB association 
(or, in the Gemini data, any object at all within our error 
ellipse) is detected in any of the optical observations — we 
therefore consider it unlikely that M74 X-l is a Galac- 
tic LMXB. Furthermore, it is unlikely that M74 X-l is a 
Galactic cataclysmic variable, since large-amplitude X- 
ray oscillations of this type are most often seen in polars 
and intermediate polars, and would then be related to 
the spin period of the white dwarf. This would result 
in highly sinusoidal periodic variation. However, the ob- 
served variability is only quasi-periodic, and in June, the 
flares are only observed toward the end of the observa- 
tion. Therefore, M74 X-l most likely resides outside our 
Galaxy. 

It is also possible that M74 X-l is a background ac- 
ti ve galactic nucl e us (A GN). Using equations 1 and 2 
of IGiacconi et all lj2001|) . and ignoring extinction from 
M74, it can be seen that zero to two sources of this flux 
would be expected on a single ACIS chip. If it were 
an AGN, it would most likely be a Narrow-Line Seyfert 
1 (NLS1): these objects are highly X-ray variable, and 
are often not radio-loud — but note that they do have 
optical counterparts, unlike M74 X-l. Taking the de- 
tection limit of my = 25 as the magnitude of the opti- 
cal counterpart, the lower limit to the X-ray-to-optica l 
flux ratio is log[F x /Fy} = 2.3 ijMaccacaro et al.l ll982). 
Here, we assume no excess absorption — which would in- 
crease this ratio — since the absorption present in the 
spectra is always consistent with the Galactic value of 
TVh- Among AGN, BL Lac objects have the highest X- 
ray to o ptical flux ratio, ra nging from \og[Fx/Fy] = 0.3 
to 1.7 ijStocke et al.lll99lj) : this is significantly lower 
than what is observed for M74 X-l. Furthermore, there 
are only a few AGN which approach the variability ex- 
hibited by M74 X-l; for example IRA S 13244-3809 
Dewangan et alJl2002tlBoller et alll997h and PHL 1092 
Brandt et al.lll999|) . Since M74 X-l's extreme variabil- 
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ity and lack of detected radio or optical counterpart all 
argue against its identification as a background AGN, we 
conclude that it almost certainly resides in M74. 

4.2. Limits on the binary system from optical data 

Assuming that M74 X-l is located in M74, the peak 
luminosity of 1.2 x 10 40 erg s" 1 implies that M74 X- 
1 most likely contains a black hole, and probably has a 
binary companion. At the distance of M74, the detection 
limit of my ~ 25 corresponds to My > —4.7; thus any 
possible comp anion must be fainter than a B OIII star 
(M v ~ -4.9. Ide Jaeer fc Nieuwenhuiizenlll987|) . 

The Gemini data with the Chandra positional error 
circle and nearby OB associations overlaid are shown 
in Figure El M74 X-l is near (but clearly not in) 
an OB association which has my = 20.6 (OB 121, see 
llvanov et al.lll"99 2fK consistent with a group of approxi- 
mately 25 O or B stars. This is interesting, since ULXs 
are often located near star-forming regi ons (see, e.g. 
iZezas fc Fahbiandl2002t iGoad et alJl2002fl. Following a 
simila r argument to that presented in IZezas fc Fabbianol 
(2002), we consider the possibility that M74 X-l is a bi- 
nary which was ejected from the nearby star cluster OB 
121. In this case, its projected separation of 2" (cor- 
responding to 85 pc) could be attained with a velocity 
in the plane of the sky of 30 km s _1 for a period of 
approximately 3 x 10 6 years. This velocity is not unex- 
pected for stellar- mass black- hole binaries, given the ob- 
served velocities of such binaries and radio pulsars within 
our Galaxy. In fact, the velocity could be significantly 
higher; the Galactic micr oquasar GRO 1655— 40 has a ve- 
locity of - 100 km s" 1 ijMirabel et al]l2002j) — adopting 
this value for M74 X-l implies that it was ejected only 
~ 8 x 10 5 years ago. The implied lifetime of ~ 10 6 
years is also compatible with the lifetime of the X-ray 
emitti ng phase of a Mrh ~ 10 black-hole X-ray bi- 
nary ijPodsiadlowski. Rappaport fc Hani 120031) . Taking 
the average luminosity of the system to be 10 erg s _1 , 
the average M is given by the equation M avo = L ave /r)C 2 . 
Since rj ~ 0.1 for black hole systems, M ~ 2 x 10~ 7 Mq 
yr _1 . Assuming that the primary is a black hole of 
10 Mq, this M implies that Mi > 8 Mq, and the optical 
data imply that Mi < 20 (the lower limit is taken 
from Figure 4 of Podsia dTowski. Rappaport fc Hanll2003l 
and the upper limit is approximately the mass of an 09V 
or BOIII star, both of which are at the luminosity limit 
for detection in the GMOS data). However, we caution 
that the peak luminosity of M74 X-l is an order of mag- 
nitude higher than the Eddington limit for a 10 M© black 
hole; this is discussed further in the following sections. 
Finally, if M74 X-l is a product of the recent star for- 
mation seen in the nearby OB associations, it is unlikely 
that it is in a system with a low-mass companion, since 
the secondary would not have had time to evolve off the 
main sequence and fill its Roche lobe. We conclude that 
M74 X-l most likely has an intermediate- or high-mass 
donor if it is related to the nearby star- forming regions. 

4.3. Interpreting the soft excess 

The spectra of M74 X-l are reasonably fit, in most 
cases, by a power-law model, although in some cases the 
addition of a soft component (such as a blackbody) is 
favored with high (> 3u) significance. Assuming that 



this component is due to an accretion disk, we consider 
the physical implications of the fitted disk blackbody, 
which has an inner disk temperature of order 0.3 keV. 
The observed inner disk temperature in the 2001 Octo- 
ber 19 "bright" data (0.29 keV) and disk blackbody lu- 
minosity (1.0 x 10 39 erg s _1 ) imply an inner disk radius 
of ~ 1.2 x 10 3 km, which, if taken to be the innermost 
stable circular orbit of a non-spinning BH, re quires a 
BH mass of ~ 140 M ijMakishima et alJl200d assum- 
ing a Schwarzschild BH.) Low-temperature disk black- 
body components have also been seen in other ULXs, 
strengthening the^ hypothesis that they contain IMBHs 

Alternatively, if M74 X-l is a stellar-mass black hole 
(which obeys the Eddington limit) then we might in- 
terpret its high luminosities as due to some sort of 
anisotropic emission process. As discussed in Section l4~Tl 
the location of the M74 X-l near OB 121 is consistent 
with an ejected black hole X-ray binary with a primary 
mass of 10 Mq. In this case, the T in of ~ 0.3 keV can 
be compared to the ~ 1 keV value found in Galactic 
BHX N when they are a ccreting at their Eddington rate 
(Makishi ma et alJ 2000). If the disk extends down to the 
last stable orbit, this lower temperature implies a lower 
mass accretion rate, T m ~ M 1//4 ijFrank. King fc Rainel 
119921 equation 5.41), which implies that M is approxi- 
mately 100 times below Eddington. Given that the ob- 
served peak disk blackbody luminosity is approximately 
the Eddington luminosity, this implies that M74 X-l is 
beamed by a factor of b ~ 0.01, where b = CI/Att. This 
is more than one would expect from simple scattering of 
the fl ux off the inner edge of a funnel shaped accretion 
disk l|King et alJl2001^ . 

It is therefore difficult to reconcile the high luminos- 
ity of M74 X-l and the low temperature of the ther- 
mal component without concluding that M74 X-l is an 
intermediate-mass black hole. However, the presence of 
a thermal component is not necessarily due to an accre- 
tion disk. Some portion of the soft flux may be due to 
another proce ss, such as X-rays s cattered in an optically- 
thick outflow ijMukai et al.l l2003). Or, it may be due to 
hot diffuse gas associated with the source. This interpre- 
tation is supported by the fact that the soft component 
is also fit well by a MEKAL model (see, e.g., Roberts et 
al. 2004 for a discussion of the physical interpretation of 
MEKAL components in ULX spectra). 

4.4. Spectral and temporal variability and the jet model 

There is some evidence that the 2001 October 19 spec- 
trum of M74 X-l hardens as it becomes brighter (see 
the discussion in Section I2.1JI . Note that similar be- 
havior has been observed in other ULXs: for exam- 
ple, a number of ULXs in NG C 4038/9 (the Anten - 
nae) harden with increasing flux llFabbiano et al][2 003L 
as do e s NGC 6946 X - ll (M F 16: IRoberts fc ColberU 
I2003|) . iFabbiano et alJ l|2003|) posit that this is due to 
an increase in the temperature of the inner disk at very 
high accretion rates. However, this model is not con- 
sistent with our observations of M74 X-l, since its spec- 
trum is consistently dominated by the power-law compo- 
nent and si nce the variability occurs on extremely short 
timescales. IRoberts fc Colbert) 1)20031) suggest two possi- 
ble models for the hardening seen in NGC 6946 X-ll — 
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that it may be due to variable emission from the base 
of a jet l|GeorganoDOulos. Aharonian fc Kirld 120021) . or 
that it may be due to ma gnetic reconnecti on events in 
the accretion disk corona ijdi Matted Il998j) . We favor 
the jet model for M74 X-l, since it also allows for the 
possibility of relativistic boosting of the X-rays and thus 
an amplification of the intrinsic variability. Relativistic 
beaming is a widely accepted explanation for many of 
the properties of face-on AGN, and several Galactic jet 
sources have also demonstrated variability that may be 
consistent with beame d X-ray emission: XTE J1650— 500 
lITomsick et all l200l and SAX J1 819.3-2525 (V4641 
Sgr. lWiinands fc van derl^l2T)0fl . 

M74 X-l exhibits flares which appear relatively equally 
spaced, and an FFT of the October observation shows 
that there is a significant broad peak in the low-frequency 
range of the power spectrum. However, there is no co- 
herent period, as can also be seen by the multiple peaks 
in Figure IT21 We therefore rule out periodic sources for 
the variability, such as an orbital or spin modulation. 
We consider the possibility of eclipses of the inner disk 
by clumpy material to be unlikely, since there is no in- 
dication of an increase in iVjj during the low-count-rate 
periods (it is consistent with the Galactic value through- 
out). It is also unlikely that the variability is due to 
jet precession, since it should then be strictly periodic, 
which it is not; furthermore, the timescales are very short 
compared with, for example, the 1 64 day precession pe- 
riod seen in SS 433 ljMargonl ll984'). However, aperiodic 
variability of this kind could be related to processes at 
the base of the jet, namely the ejection of relativistic 
plasma into the jet. While such processes are not well 
understood, there is an observed correlation between ac- 
tivity in the inner disk and at the base of the jet in the 
Galactic microquasar GRS 1915+105, which happens on 
timescales similar to what is observed in M74 X-l (see 
Section 1431 for further discussion). 

In each of its three observations, M74 X-l is extremely 
variable and exhibits powerful X-ray flares — its luminos- 
ity fluctuates by an order of magnitude on timescales of 
half an hour. This is suggestive of beaming due to a rel- 
ativistic jet: as seen in BL Lac objects (commonly called 
blazars), relativistic beaming along the line of sight leads 
to rapid, high-amplitude variability. M74 X-l may, in 
fact, be a "microblazar" — a microquasar whose relativis- 
tic jet is beaming the X-ray emission towards us along 
the line of sight. 

4.5. Comparison with other ULXs and Galactic 
microquasars 

Of Galactic sources, M74 X-l closely resembles micro- 
quasars in both its temporal variability and energetics. A 
link between microqua sars and ULXs has been posited by 
a number of authors |K ordin^ 1 Falcke fc Markofj 120021 
IGeoreanopoulos. Aharonian fc Kirkll2002l for example'! ; 
and recently, r adio emission was d iscovered from a ULX 
in NGC 5408 ijKaaret et al.ll2003|) . This radio emission 
is almost certainly due to a relativistic jet, and given the 
high radio flux, it is likely beamed, with a Doppler fac- 
tor S = 7 _1 (1 — ficosO)" 1 > 5.8. The authors point out 
that this is consistent with the known ejecta velocities 
of GRS 1915+105, which would produce such a Doppler 
factor given a jet axis tilted up to 10° with respect to the 
line of sight (they further note that approximately 1 in 70 



extragalactic microquasars should be thus aligned) . The 
X-ray spectrum of the NGC 5408 source is also consis- 
tent with beamed emission, as it is well fit with a broken 
power-law model (but is also reasonably fit by a disk 
blackbody of 0.1 keV plus a power-law). 

GRS 1915+105 is very similar to M74 X-l in its tempo- 
ral characteristics as well as its persistence. It has been 
active since its discov ery in 1992, exhibiting a wide range 
of behavior (see, e.g.. IMuno. Morgan fc Remillardlll999t 
Mar kwardt. Swank fc T aam 1999). In some observa- 
tions, GRS 1915+105 has b een observed to vary from 5 x 
1(F to 3.5 x 10 39 erg s" 1 ijMarkwardt. Swank fc Taanl 
119991 assuming D = 12.5 kpc). This can be compared to 
M74 X-l, which has been active since 1980, and varies 
from - 5 x 10 38 to 1.2 x 10 40 erg s -1 . Although the low 
disk temperature and higher peak luminosities observed 
in M74 X-l favor a larger black hole mass, the flares seen 
in M74 X-l occur on the same timescales as some flar- 
ing behavior seen in GRS 1915+105, suggesting a similar 
dynamical timescale and thus similar physical scale. 

Although much of the 0.3-10.0 keV X-ray emis- 
sion from GRS 1915+105 is thought to come from its 
disk, there is also evidence for a correlation between 
disk emission and jet formation ((Mirabel et al.l Il998t 
lEikenberrv et al1 fT998) . Therefore, the timescale for the 
variability at the base of the jet is likely similar to what 
is seen in X-ray observations. The jets of GRS 1915+105 
are not directed toward our line of sig ht, but are inclined 
by 70° l|Mirabel fc Rodriguez! IT9941. If a jet were di- 
rected toward us, we would expect to see relativistically 
boosted emission varying on similar timescales to what is 
currently observed: we propose this scenario as a possi- 
ble explanation for the high-amplitude flares seen in M74 
X-l. 

Another prominent feature of Galactic microquasars is 
their activity at radio wavelengths, which includes vari- 
ability on timescales of tens of minutes to several hours. 
The radio flux of GRS 1915+105 has been observed to 
be as high as 430 mJy at 3.6 cm l)Rodriguez fc Mirabel 
11997ft . but at the distance of M74 this would only be 0.8 
/iJy — unfortunately, not detectable by any current radio 
telescopes. Even if M74 X-l were as radio-bright as the 
source in NGC 5408, at the distance of M74, it would 
be at the limit of detectability. As mentioned above, a 
search for radio emission in a 1.425 GHz VLA observa- 
tion did not find emission at the location of M74 X-l 
with an upper limit of 80 /iJy. 

5. SUMMARY AND CONCLUSIONS 

M74 X-l is a highly variable X-ray source whose 0.3- 
10.0 keV luminosity is seen to vary from approximately 
5 x 10 38 to 1.2 x 10 40 erg s _1 in several thousand sec- 
onds. It was detected by Einstein in 1980, and observed 
twice with Chandra and once with XMM-Newton. The 
latter three observations of M74 X-l were each separated 
by approximately four months. There have also been a 
number of optical observations and a radio observation 
of the region, but as of yet, no counterpart has been 
found. Since M74 X-l has no optical counterpart and 
is extremely variable, it is most likely in M74 (and thus 
a ULX), and may be related to the nearby OB associa- 
tion. In addition, its proximity to an OB association (a 
relatively common occurrence for ULXs) implies that it 
is likely a high- or intermediate-mass system. 



10 



The spectrum of M74 X-l is consistent with a low tem- 
perature (~ 0.3 keV) disk blackbody component which 
has been observed in other ULXs. This component is 
often interpreted as arising from an accretion disk; its 
low temperature then implies an inner disk radius of 
1.2 x 10 3 km and therefore a black hole primary of mass 
140 Mq. However, in many ways M74 X-l resembles 
stellar-mass Galactic microquasars. In this scenario, the 
very high X-ray luminosity and high-amplitude rapid 
variability would be produced by beaming from a jet 
which is aligned close to our line of sight. If M74 X-l 
is, in fact, a stellar-mass black hole of mass ~ 10 Mq, it 
becomes difficult to explain the low temperature compo- 
nent in terms of an accretion disk, and it may be from a 
different physical source. 

M74 X-l is an intriguing object. It appears to be a 
black hole intermediate- or high-mass X-ray binary where 
the mass of the black hole remains an open question. 
We favor an interpretation in which M74 X-l harbors 
a stellar mass black hole and the high observed X-ray 
luminosity is produced by beaming. However, we cannot 



rule out an intermediate-mass black hole interpretation 
with the current data. 
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